Diffractive structure function from fits with higher twist 
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Abstract 

We make predictions for the diffractive longitudinal structure function 
FjP to be measured at HERA, based on DGLAP based fits of diffrac- 
tive parton distributions with twist-4 contribution. This contribution 
describes diffractive qq production from longitudinal photons and sig- 
nificantly changes predictions for F£ obtained in pure DGLAP analy- 
ses. 

1 Introduction 

We are interested in diffractive deep inelastic scattering (DDIS) at HERA which provide a very 
interesting example of processes with a clear interplay between hard and soft aspects of QCD 
interactions. In these processes, a diffractive system is formed which is separated in rapidity from 
the scattered proton. The most important observation made at HERA is that diffractive processes 
in DIS are not rare and constitute up to 15% of all deep inelastic scattering events [1-4]. 

After integration over the proton azimuthal angle, the diffractive cross section is given in 
terms two structure functions, and , which depend on four variables: Bjorken-x, photon 
virtuality Q 2 and two additional variables 

Q 2 + M 2 -t 2 
X]P= Q2 + w 2 ' t = (p-P) • (!) 

Here M is mass of the diffractive system, W is invariant energy of the gamma-proton system 
and p, p' are incident and scattered proton momenta. In our analysis, the diffractive structure 
functions are decomposed into the leading twist-2 and higher twist contributions 

F 2 D L = F^ w2) + F^ tw4) + ... . (2) 

In the Bjorken limit, the leading part depends logarithmically on Q 2 while the twists part is sup- 
pressed by additional power of l/Q 2 . However, this contribution dominates over the twist-2 one 
for small diffractive masses, M 2 <C Q 2 , playing especially important role in DDIS. Physically, 
the twist-4 contribution is given by diffractively produced qq pairs from longitudinally polar- 
ized virtual photons. The effect of this contribution is particularly important for the diffractive 
longitudinal structure function F® which is supposed to be determined from the HERA data. 
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2 Twist-2 contribution 

The twist-2 contribution is given in terms of diffractive partem distributions (DPD) through stan- 
dard collinear factorization formulae [5-8]. In the next-to-leading logarithmic approximation we 
have 

F? (tw2 \x, Q 2 ,x F , t) = Sd + ^ {C! ®S d + C 2 g ® G D ) (3) 
F^ tw2 \x,Q 2 , XlP ,t) = g{C7f ®S D + C% ®G D ) (4) 

S G 

where a s is the strong coupling constant and C 2 ' L are coefficients functions, known from inclu- 
sive DIS [9, 10]. The integral convolution is performed for the longitudinal momentum fraction 
and reads 

(C®F)(J3) = [ 1 dzC(P/z)F(z). (5) 

Notice that in the leading order, when terms proportional to a s are neg lected, F? {tw2) = 0. The 
functions So and Go are built from diffractive quark (q* D ) and gluon {go) distributions 

N f 

S D (x,Q 2 , XlP ,t) = J^e} (3 {q f D ((3,Q 2 ,x 1P ,t) +q f D (P,Q 2 ,xjp,t)} (6) 
/=i 

and 

G D (x,Q 2 ,x P ,t) = (3g D (l3,Q 2 ,x 1P ,t). (7) 
From these equations we see that the new variable 

(3 = ^ = g^W 2 (8) 

plays the role of the Bjorken variable in DDIS. 

The DPD evolve with Q 2 with the DGLAP evolution equations [1 1] for which (xjp, t) are 
external parameters. In this analysis we assume Regge factorization for these variables: 

q f D (P, Q 2 ,xjp, t) = f P (x F , t) q f jp(P, Q 2 ) (9) 
g D {f3,Q 2 ,X]p,t) = f P (x P ,t) g F (/3,Q 2 ). (10) 

The motivation for such a factorization is a model of diffractive interactions with the pomeron 
exchange [12]. In this model fjp is the pomeron flux 

fp{xp,t) = N ^ 2 X]p w , (11) 
where ajp(t) = ajp(0) + a'jpt is the pomeron Regge trajectory and 

Fl{t) = i£(0)e- B "l*l, (12) 



Fig. 1: Twist-4 contribution from longitudinally polarized photons. Two gluons here model the pomeron exchange 
which is later unitarized and effectively described by the dipole cross section, see the text below. 



is elastic formfactor which describes the pomeron-proton coupling. The diffractive slope is taken 
form HERA data, B D = 5.5 GeV~ 2 , and Ff>(0) = 54.4 GeV~ 2 [7]. For the pomeron trajectory, 
a'jp = 0.25 GeV 2 but the intercept ajp(0) is fitted to the analyzed data. The diffractive quark 
distributions are flavour independent, thus 

<£(A Q 2 ) = ¥ P (P, Q 2 ) = Q 2 ) d3) 

where Nf is a number of active flavours and Sjp is singlet distribution. We fit £_p(/3, Q 2 ) and 
Gjp(/3, Q 2 ) at an initial scale Ql = 1.5 GeV 2 to diffractive data from HERA, using the DGLAP 
evolution equations in the next-to-leading order approximation. We also include charm quark 
contribution into the analysis [13]. 



3 Twist-4 contribution 

The twists contribution describes diffractive production of the qq pairs from longitudinally po- 
larized virtual photons, see Fig. [T] Although formally suppressed by l/Q 2 , this contribution 
dominates over the twist-2 contribution for small diffractive masses, M 2 <C Q 2 (or [5 — > 1) 
[14-16]. This is why it cannot be neglected in the analysis of diffractive DIS data. 

We used the following form of the twist-^- contribution which has to be added to the 
diffractive structure functions p^ tw2 ^ an( j pD(tw2) 



4 3 

F r D(tw4) = - e~ BD ^ V e 2 ( f dk 2 k2/Q2 - <b 2 (14) 

V 1 i-/?^ 2 



with 




o 



kr Jo(kr) a(xjp, r) 



(15) 



where Kq and Jo are Bessel functions. Strictly speaking, formula (fT4l contains all powers of 
l/Q 2 but the twist-4 part, proportional to l/Q 2 , dominates. The function a(xp,r) in eq. (fT5T ) 
is the dipole-proton cross section which describes the interaction of a quark-antiquark dipole of 
transverse size r with the proton. Following [18], we choose 



where Q 2 (xp) = {xp/xq)~ x GeV 2 is a saturation scale which provides energy dependence of 
the twist-^1 contribution. The dipole cross section parameters, uq = 29 mb, xq = 4 • 10~ 5 and 
A = 0.28, are taken from [18] (Fit 2 with charm). This form of the dipole cross section provides 
successful description of the inclusive and diffractive data from HERA. 

In addition to the twist terms, we also consider a reggeon contribution, described in detail 
in [13], which improves fit quality through better dependence on xjp. 

4 Fit results 

In our analysis, we use diffractive data on (or reduced cross section a®) from the HI [3,4] 
and ZEUS [1,2] collaborations. These data were obtained in different kinematical regions, using 
different methods of their analysis, thus we decided to analyse them separately. For each data set 
we performed two fits: with the twists present and without this term (pure DGLAP fits). In this 
way, we obtained two sets of diffractive parton distributions which allow us to make predictions 
for the longitudinal structure function FjP. A full discussion of fit details is given in [13]. 

The diffractive parton distributions (DPD) from the analysis of HI data are shown in Fig. [2] 
They are given in terms of the pomeron parton distributions which can be multiplied by the 
pomeron flux /(xjp, t) to obtain the DPD. We see that the singlet quark distributions from the 
two discussed fits are practically the same, while the gluon distributions are different. The gluon 
from the fit with twist-4 is stronger peaked near ss 1 than in the twist-2 fit. This somewhat 
surprising result can be understood by analyzing the logarithmic slope dF® /dlnQ 2 for fixed f3. 
In the leading logarithmic approximation, we have from the DGLAP equations 



where the negative term sums virtual corrections. For large (3, the measured slope is negative 
which means that the negative term must dominate over the positive terms describing real emis- 
sions. The addition of the twist-4 contribution to F®, proportional to l/Q 2 , contributes negative 
value to the slope which has to be compensated by a larger gluon distribution near j3 w 1 in order 
to describe the same data. 

In Fig.|3l we show the diffractive structure functions resulting from the determined parton 
distributions. As expected, F 2 D is practically the same in both fits. However, the F® curves are 
significantly different due to the twist-^1 contribution (shown as the dotted lines). Let us em- 
phasised that both sets of curves were found in the fits which well describe the existing data, 



&(x]p,r) = a {1 - exp {-r 2 Q 2 s {x p))} 



(16) 




(17) 
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Fig. 2: Pomeron partem distributions: quark singlet (left) and gluon (right) from the HI data. Solid lines: twist-2 fit; 
dashed lines: fit with twist-4 term. 



especially in the region of large (3 where twist-4 is important. Thus, an independent measure- 
ment of FjP in this region would be very important for confirmation of the QCD approach to 
diffraction. 

We summarize the effect of the twist-4 contribution in Fig.|H showing the predictions for 
the diffractive longitudinal structure function. Ignoring this contribution, we find the two solid 
curves coming from the pure twist-2 analyses of the HI (upper) and ZEUS (lower) data. With 
twist-4, the dashed curves are found, the upper curve from the HI data and the lower one from 
the ZEUS data analyses. There is a significant difference between these predictions in the region 
of large 0. We believe that this effect will be confirmed by the forthcoming analysis of the HERA 
data. 



DSF (H1) 




Fig. 3: Diffractive structure functions F^ (left) and (right) from fits to the HI data for xp — 10 3 . Solid lines: 
twist-2 fit; dashed lines: twist-(2+4) fit; dotted lines: twist-4 contribution. 



5 Summary 

We performed the analysis of the diffractive data from HERA, determining diffractive parton 
distributions. In addition to the standard twist-2 formulae, we also considered the twists contri- 
bution. This contribution comes from the qq diffractive production from longitudinally polarized 
virtual photons and dominates for M 2 <C Q 2 (large j3). The twist-4 contribution leads to the 
diffractive gluon distribution which is stronger peaked at j3 « 1 than the gluon distribution from 
the pure twist-2 fits. 

The main result of our analysis is a new prediction for the longitudinal diffractive structure 
function . The twists term significantly enhances in the region of (3 > 0.6. A measure- 
ment of this function at HERA in the region of large f3 should confirm the presented expectations 
which are based on perturbative QCD calculations. 
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Fig. 4: Predictions for F° (3> for x F = 10~ 3 and Q 2 = 10 GeV 2 from fits with twist-4 to the HI (upper dashed 
line) and ZEUS (lower dashed line) data. The solid lines show predictions from pure twist-2 fits to the same data: 
from HI (upper) and ZEUS (lower). 
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